arXiv: 1501.06608vl [cond-mat.mes-hall] 26 Jan 2015 


Coulomb bound states and resonances due to groups of Ca dimers adsorbed on 

suspended graphene 


Alireza Saffarzadeh^’and George Kirczenow^ 

^Department of Physies, Payame Noor University, P.O. Box 19395-3697 Tehran, Iran 
^Department of Physies, Simon Fraser University, Burnaby, British Columbia, Canada V5A 1S6 

(Dated: Wednesday 28^^ January, 2015) 

The electronic bound states and resonances in the vicinity of the Dirac point energy due to the 
adsorption of calcium dimers on a suspended graphene monolayer are explored theoretically using 
density functional theory (DFT) and an improved extended Hiickel model that includes electrostatic 
potentials. The Mulliken atomic charges and the electrostatic potentials are obtained from DFT 
calculations and reveal charge transfer from the Ca dimers to the graphene which is responsible 
for the emergence of resonant states in the electronic spectrum. The number of resonant states 
increases as the number of adsorbed dimers is increased. We hnd a bound “atomic-collapse” state 
in the graphene local density of states, as has been observed experimentally [Wang et ah, Science 
340, 734 (2013)]. We hnd the formation of the atomic-collapse state and its population with 
electrons to require fewer adsorbed Ca dimers than in the experiment, possibly due to the different 
spacing between dimers and the dielectric screening by a boron nitride substrate in the experiment. 
We also predict the onset of Filing of a second atomic-collapse state with electrons when six Ca 
dimers are adsorbed on the suspended graphene monolayer. Experiments testing these predictions 
would be of interest. 


I. INTRODUCTION 

The electronic properties of graphene have attracted 
a great deal of research due to their fundamental inter¬ 
est and promising applications [T] which are related to 
graphene’s linear dispersion near the Dirac point energy 
[21 Hj. In this regard, the problem of charged impurities 
in graphene has recently received considerable attention 
[HE]. These impurities, as well as adsorbates, can in¬ 
duce resonant states in the vicinity of the Dirac point 
energy that play an important role in a number of prop¬ 
erties of graphene nanostructures, such as transport 0- 
HS], chemical sensing inHis], controlled doping [IZIIIHI 
and magnetism [T9]-[29]. 

It was predicted theoretically that when the Coulomb 
potential strength for heavily charged impurities in 
graphene exceeds a certain critical value, an infinite fam¬ 
ily of Rydberg-like quasi-bound states appears abruptly 
[6]. Therefore, as an atomic collapse (fall to the cen¬ 
ter effect) was predicted for isolated atoms with highly 
charged nuclei, one can expect a similar situation with 
highly charged impurities in graphene [6]. Such a phe¬ 
nomenon was observed recently in the scanning tunnel¬ 
ing microscopy (STM) experiments of Wang et al. [30] 
in ultra-high vacuum. Using atomic manipulation, they 
assembled groups of dimers of calcium ions on the sur¬ 
face of graphene placed on a boron nitride (BN) sub¬ 
strate. Ca dimers were added to the group, one by one, 
until the total charge transferred from the Ca to the 
graphene surpassed a critical threshold m- The mea¬ 
sured spatial and energetic characteristics of the elec- 
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tronic states around the Ca-dimer cluster showed the 
emergence of an oscillation in the local density of states 
(LDOS) that manifested itself as a resonance above the 
graphene Dirac point as the number of dimers in the clus¬ 
ter was increased m- When the resonance shifted below 
the Dirac point, the quasi-bound state was interpreted 
as the atomic-collapse eigenstate m- In support of this 
interpretation of the experimental results, calculations 
were performed m to simulate the LDOS of graphene 
near charged impurities by means of a parameterized 
model for charged impurities and the two-dimensional 
(2D) continuum Dirac model HHS). An ab initio density 
functional theory (DFT) calculation was performed [30| 
for a single Ca dimer adsorbed on graphene to determine 
a reasonable value for the model fitting parameter. These 
model calculations m were able to reproduce the main 
effect observed in the experiment [30] , namely, formation 
of the the bound electronic state that is the analog of 
the atomic-collapse phenomenon. However, a more de¬ 
tailed theoretical treatment of the graphene electronic 
states associated with groups of Ca dimers adsorbed on 
the graphene is clearly desirable. Specifically, ab initio 
DFT calculations of the charge transfer and electrostatic 
potentials associated with more than just one adsorbed 
Ca dimer (the case treated with DFT in Ref. [3Q|) . and 
calculations of the resonant and bound electronic states 
based on atomistic graphene models, as distinct from the 
continuum model considered in Ref. EOl have not been 
available. Furthermore, the previous theoretical study 
m did not address the dependence of the spatial and 
energetic characteristics of the resonant states on the spa¬ 
tial arrangement of the Ca dimers on graphene. 

In this paper, we explore the resonance features seen 
in the STM experiment m by means of DFT calcula¬ 
tions of the relaxed geometries, charge transfer between 
the Ca dimers and graphene, and the associated electro- 
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static potentials for groups of 1-6 Ca dimers adsorbed 
on suspended graphene. In order to model the electronic 
structures for groups of Ca dimers adsorbed on graphene 
substrates with large numbers of carbon atoms atomisti- 
cally, we use the (tight-binding) extended Hiickel model 
of quantum chemistry [3TH33] augmented by the inclu¬ 
sion of electrostatic potentials obtained from DFT calcu¬ 
lations. The LDOS for graphene in the presence of 1-6 
Ca dimers is calculated and the emergence and evolu¬ 
tion of resonances and their spatial characteristics with 
increasing the numbers of Ca dimers are discussed. In 
addition, we compare the Dirac point resonances for six- 
dimer clusters (not addressed in the STM experiment of 
Wang et al ED]) with two different dimer arrangements. 
Our calculations show that the resonance features in the 
LDOS are qualitatively consistent with the experiment 
and that the emergence of these resonances is due to the 
charge transfer from Ca atoms to graphene monolayer 
and not due to the atomic orbitals localized on the Ca 
atoms. Furthermore, we find that the arrangement of the 
Ca dimers on the graphene can affect the energies and os¬ 
cillations of the resonances somewhat at short distances 
from the centers of the Ca-dimer clusters. 

However, we predict formation of the atomic-collapse 
state and its population with electrons to require fewer 
Ca dimers in the cluster adsorbed on the graphene than 
was required in the experiment of Wang et al. [30] . We 
attribute this difference to dielectric screening by the BN 
substrate underlying the graphene in the experiment and 
to the larger spacing between the Ca dimers in the clus¬ 
ters studied in the experiment. m Furthermore, for 
six adsorbed dimers (a case not addressed in Ref. |30|), 
we predict a second atomic-collapse state to also become 
partly populated with electrons. Experimental tests of 
these predictions for compact clusters of Ca dimers on 
suspended graphene would be of interest. 

The paper is organized as follows. In Sec. II, we 
present the model Hamiltonian for Ca dimers adsorbed 
on graphene, derived from extended Hiickel theory aug¬ 
mented by the inclusion of electrostatic potentials ob¬ 
tained from DFT calculations. We also discuss the pro¬ 
cedure for designing large Ca dimer/graphene clusters for 
the purpose of eliminating edge effects in the electronic 
states close to the graphene Dirac point. The Mulliken 
atomic charges and the electrostatic potentials on the 
graphene monolayer, obtained from DFT calculations in 
the presence of Ca dimers, are discussed in Sec. HI.A. 
The calculated LDOS within an energy window close to 
the Dirac point of graphene at different distances from 
the centers of the Ca-dimer clusters and discussions of 
the resonant states, the atomic-collapse state, and state 
filling with electrons are presented in Sec. III.B. In Sec. 

III. C, the dependence of resonance states and their filling 
with electrons on the spatial arrangement of the dimers 
is examined by considering a six-dimer cluster and a six- 
dimer ring on the graphene monolayer. Finally, in Sec. 

IV, we conclude with a general discussion of our results. 



FIG. I: (Color online) (a) The optimized geometrical struc¬ 
ture of a single Ca dimer on a graphene disk of 150 carbon 
atoms passivated by 30 hydrogen atoms, (b) Side view of (a), 
but only with 54 C atoms which shows the (< 0.31 A) out-of- 
plane (downward) distortion for the carbon atoms beneath the 
Ca dimer. For the sake of clarity, the C-C bonds are not shown 
here, (c) Two, (d) three, (e) four, and (f) five Ca dimers on a 
graphene ribbon extended along the ^-direction and designed 
based on the optimized geometry shown in (a), and sand¬ 
wiched between two semi-infinite leads extended along the 
x-direction. Only a segment of each central ribbon, referred 
to as the “dimer-graphene cluster,” is shown here. 


II. THEORY 

We carried out ab initio geometry relaxations based on 
DFT for a single calcium dimer on a honeycomb graphene 
lattice using the GAUSSIAN 09 software package with 
the PBEhIPBE hybrid function and the 6-3IG** basis 
set m- Eor this purpose we adopted a graphene disk 
consisting of 150 carbon atoms passivated at the edges 
with 30 hydrogen atoms, the Ca dimer being bonded to 
the graphene disk at its center. Although a finite-size 
graphene cluster was used for the optimization calcula¬ 
tions, the bond lengths between neighboring C atoms 
were well converged with increasing distance from the 
Ca, which indicates that the number of 150 C atoms is 
large enough to model the graphene monolayer. Since 
all of the atoms were allowed to relax freely, we found a 
maximum in-plane distortion of 0.02 A between carbon 
atoms and a maximum out-of-plane distortion of 0.31 A 
(away from the Ca) due to the presence of Ca dimer on 
the graphene. The final optimized geometry obtained in 
this study is shown in Eig. 1(a). The distance between 
the two Ca atoms of the dimer is 3.80 A and the center 
of the dimer is located at 2.4 A above the plane defined 
by the C atoms far from the Ca which is in agreement 
with the recent ab initio DET calculations in a gener¬ 
alized gradient approximation for a Ca dimer adsorbed 
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on graphene in supercells with periodic boundary condi¬ 
tions m- The positions of the Ca atoms in the graphene 
plane (relative to the center of the graphene disk) are 
Ri = —O.SSAx + l.bSAy^ and R 2 = — 0.20Ax — 2.2lA^, 
indicating that the Ca atoms are in nonequivalent posi¬ 
tions. 

Since the adsorbed Ca dimers on graphene behave as 
charged impurities, the electrostatic properties such as 
spatial charge distributions and electrostatic potentials 
at each atomic site should be computed by means of 
DFT calculations. For this reason, we designed a larger 
graphene hexagonal disk similar to that shown in Fig. 
1(a) but with 1370 C atoms (not shown here). Then, 
we transferred the atomic coordinates of the relaxed Ca 
dimer along with those of the 54 closest C atoms sur¬ 
rounding the dimer and belonging to the optimized clus¬ 
ter in Fig. 1(a) to the locations of each of the Ca dimers 
of the larger disk to produce the structures with multiple 
Ca dimers shown in Fig.l(c)-(f). This procedure is justi¬ 
fied since the distortion of the graphene lattice induced 
by an adsorbed Ca dimer has a shorter range than the 
distances between the Ca dimers in the clusters consid¬ 
ered in this work (see Fig. 1(b)).The cluster geometries 
formed in this way were used to find the Mulliken atomic 
charges and the electric potentials on the carbon atoms 
of the graphene and the calcium atoms by means of DFT 
calculations. 

The main difficulty in recognizing the physical reso¬ 
nant states induced by Ca dimers, especially at large dis¬ 
tances from the dimer region where the influence of res¬ 
onant states is not strong, is the mixing of these states 
with fluctuations in the electronic states around Dirac 
point energy due to edge effects in the model system. In 
order to suppress these quantum size effects as much as 
possible, a large rectangular graphene sheet (called here 
the “central ribbon”) consisting of 3940 carbon atoms 
with length 41.2 A in the x-direction and width 240 A in 
the ^-direction was designed. Then all the relevant opti¬ 
mized coordinates and the atomic electrostatic potentials 
obtained from the DFT calculations for the 1370 C atom 
disk with adsorbed Ca dimers were applied to the appro¬ 
priate atoms in this ribbon. Next, we attached two semi¬ 
infinite graphene ribbons (“leads” hereafter) with width 
240 A to the left and right sides of the central ribbon 
to produce an infinite ribbon (in the x-direction) with 
width 240 A. It is assumed that, the electrostatic poten¬ 
tial values for all carbon atoms located far away from the 
dimers region converge to a fixed value. Therefore, for 
the remaining carbon atoms in the central ribbon and 
for all the atoms in the semi-infinite leads a cut-off in 
the electrostatic potential was used. A discussion of this 
point will be given in the next section. 

According to the description given above, the Hamil¬ 
tonian of the system which is partitioned in three blocks, 
i.e. the central ribbon C and the two semi-infinite leads 
at its left (L) and right {R) sides, can be written as 

( 1 ) 


where the Hamiltonians of the central ribbon, Hc^ and 
the leads, Hl,r^ are described by a tight-binding model 
derived from the extended Hiickel theory in a basis of C 
25 and 2p, and Ca 45 and 4p valence orbitals, that can 
be written as 

^iE{C,L,R} ^ ^ "t" '7i,a/3 T H.C.) , 

cx a,(3 

( 2 ) 

where dl^{dia) is the creation (annihilation) operator for 
an electron in the ath atomic valence orbital ^pio^. Within 
extended Hiickel theory eia is the experimentally deter¬ 
mined valence orbital ionization energy, and 7i,a/3 is the 
off-diagonal matrix element between valence orbitals pJia 
and In Eq. (IT|), VL(i?)c describes the coupling ma¬ 
trix between the lead L (R) and the central ribbon C. 
Since the basis set used in extended Hiickel theory is 
nonorthogonal, the orbital overlap can 

be non-zero and will be included in this study [35l |36] . 
To include the atomic electrostatic potentials obtained 
from our DFT calculations in the present model, the di¬ 
agonal and off-diagonal matrix elements in the Hamil¬ 
tonians are modified as ^ + Ui{xa^ya) and 

li,aP -X Xi,al3 + \[Ui{Xa,ya) + 2//3)]S'i,a/3 where 

XaiVa) is the X (y) component of the position of atom 
with orbital a, and U{x^^ya) is the electric potential en¬ 
ergy (EPE) for an atom at point R = {xa,ya)- Note 
that the EPE values that will be discussed below, are 
the same for all valence orbitals localized on each given 
atom of the dimer-graphene system. 

The Green’s function of the central ribbon coupled to 
the graphene leads can be written as [371138] 

Gc(6) = -He- ^Ue) - E^(e)]-1 , (3) 

where, z = e ^ if] with y being a positive infinitesimal 
number, and T^{l,r} ai'c the self-energy terms due to the 
semi-infinite leads, each of which is considered as a semi¬ 
infinite stack of principal layers (ribbons) with nearest- 
neighbor interactions [39] and can be defined as [40| 

^LiR){^) = , ( 4 ) 

where = yLiR)c-^SL(R)c with Sl(r)c being the 

overlap matrix between lead L{R) and the central ribbon 
G. gLiR) = {zSl(r) - 77L(i?))“^ represents the surface 
Green’s function of the lead L{R) that can be written in 
terms of appropriate transfer matrices and computed us¬ 
ing the iterative procedure introduced by Lopez-Sancho 
et al. [381440] . Accordingly, we can calculate the local 
density of states (LDOS) for the nth site in the central 
ribbon from the diagonal elements of Gc as 

pn(e) = --Im[Gc(e)S'c]n,n . (5) 

TT 

Note that in all LDOS graphs that will be presented 
in this work, the average of LDOS associated with two 
neighboring carbon atoms in a single graphene unit cell 
is taken. 


H = Hc^Hl^Hr^Vlc + Vrc 
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FIG. 2: (Color online) Distribution of Mulliken atomic 
charges (in units of |e|) obtained from DFT for (a) one dimer 
and (b) five dimers on graphene monolayer. The symbol 0 
shows the location of Ca atoms in each dimer. The Ca atoms 
are positively charged after charge transfer from the dimers 
to a graphene monolayer. 



FIG. 3: (Color online) Distribution of atomic EPE for (a) one 
dimer and (b) five dimers on a graphene monolayer. The black 
dots in (a) and (b) show the locations of the carbon atom 
nuclei, (c) The averaged EPE, U{R), as a function of distance 
from the center of dimer-graphene clusters composed of 1-6 
Ca dimers. Note that there are two clusters with different 
spatial arrangements of six Ca dimers, as discussed in Sec. 

me. 


III. RESULTS AND DISCUSSION 

A. Mulliken charges and electrostatic potentials 

Adsorption of Ca dimers on the graphene plane results 
in a charge transfer between the Ca atoms of the dimers 
and the C atoms of the graphene. To explore the nature 
of this charge transfer and the polarity of charges, the 
Mulliken atomic charges were computed by DFT for all 
of the dimer-graphene clusters. The results for the case 
of one and five dimers are shown in Figs. 2(a) and 2(b), 
respectively. Figure 2 shows that electronic charge is 
transferred from the Ca dimers into graphene states and 
the calcium atoms become positively charged, in agree¬ 
ment with the experiment [30] . In addition, the negative 
localized charges, shown by light green circles around the 
single dimer in Fig. 2(a), form a stadium-shaped region 
similar to the charge density seen in the supplementary 
materials in Ref. [30]. From the charge distributions, it 
is obvious that the Mulliken atomic charges of the car¬ 
bon atoms in the neighborhood of each single calcium 
atom are typically more negative than those of the car¬ 
bon atoms far away from the dimers, screening the bare 
charges on calcium atoms. Some evidence of Friedel-like 
charge oscillations in the graphene is also visible, espe¬ 
cially in Fig. 2(a). Importantly, increasing the number of 
dimers on the graphene sheet does not appreciably affect 
the negative-charge accumulations on the carbon atoms 
located in close vicinity to the Ca atoms, while in the 
case of a five-dimer structure, shown in Fig. 2(b), the 
concentration of negative charges on the carbon atoms 
away from the dimers is higher than for the one-dimer 
case in Fig. 2(a). The average charge transferred per Ca 
atom to the graphene depends on the number of dimers 
adsorbed on the graphene. From the calculated values 
of the Mulliken charges we found that as the number 
of dimers in a cluster was increased from one to five, 
the average charge transferred to the graphene per Ca 
atom in the cluster changes as -0.83|e|, -0.72|e|, -0.70|e|, 
-0.69|e|, and -0.69|e|, respectively. That is, somewhat 
less than one electron per Ca atom is transferred to the 
graphene and the amount of charge transferred per Ca 
atom decreases as the number of Ca atoms in the cluster 
increases, as may be expected due to the repulsive na¬ 
ture of the interaction between the electrons transferred 
to the graphene. 

The spatial distributions of the EPE, U{x^y) = 
— |e|AU, for two typical Ca-dimer clusters composed of 
one and five dimers are shown in Figs. 3(a) and 3(b). 
Here, e is the charge of electron and AV = Vc{x^y) — Vq 
is the difference between electric potential Vc for an atom 
at point R in the dimer-graphene cluster containing 1370 
carbon atoms, and Vq due to an atom of the same type 
as that at point R in the cluster but in free space. R is 
measured with respect to the center of each of the dimer- 
graphene clusters shown in Fig. 1. All of the electric 
potentials obtained from DFT were computed with the 
GAUSSIAN 09 package. Note that only the EPE at the 
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sites of carbon atoms of the graphene are shown in Fig. 
3. The EPE at the location of carbon atoms clearly re¬ 
veals the influence of adsorbed dimers on the graphene in 
both one- and five-dimer clusters. The atomic sites with 
negative (positive) EPE values represent positive (nega¬ 
tive) charges, which reflects the charge transfer from Ca 
dimers to carbon atoms of the graphene, as expected. Eor 
R > 15A away from the dimers region in the x direction, 
the EPE shows a convergence behavior due to the pres¬ 
ence of screened charges on the Ca dimers. In order to get 
a better understanding of the behavior of EPE as a func¬ 
tion oi R= {x‘^ + ^^) 2, we computed the average of EPE 
within each Ai? = lA annulus for the dimer-graphene 
systems composed of 1-6 Ca dimers (a discussion for a 
six-dimer cluster and a six-dimer ring will be given in 
Section IIIC). The results are shown in Eig. 3(c). The 
averaged EPE, U{R)^ for a single Ca dimer on graphene 
rises monotonically with distance R from the cluster cen¬ 
ter, as expected for the screened Coulomb potential due 
to a positively charged dimer. As the number of Ca 
dimers increases from one to five the electric potential 
in the interval 5A< R <25A becomes deeper due to the 
charge increase in the region of the dimers, which induces 
more quasi-bound states and stronger resonances in the 
electronic spectra of the dimer-graphene system. It is 
clear that the asymptotic behavior of U(R) is almost lin¬ 
ear for small and large distances within the range shown 
in Eig. 3(c). Eor example, in the case of a cluster with 
a single Ca dimer, the expressions U{R) = O.ISR — 0.45 
and U{R) = 0.025i? + 0.76 {U is in units of eV and R is 
in A) can be fitted to the numerical results in the range 
of R < 5Aand 15A<R<25Arespectively. 

The finite size of our 1370-carbon atom cluster used in 
the electrostatic potential calculations results in potential 
fluctuations at the edges of the cluster due to the presence 
of electronic states located at the edges. Accordingly, in 
the calculations reported below, we introduce a cutoff in 
the potential profile to lessen the influence of such edge 
effects as much as possible. The EPE cutoff value, chosen 
in all calculations that follow, is 1.2 eV. 


B. Resonances in 1-5 Ca dimer clusters 

The inclusion of the electric potential in our calcula¬ 
tions has crucial effects on the calcium-dimer-induced 
Dirac point resonant states of the graphene monolayer 
that will be studied here. Before presenting the main re¬ 
sults for the LDOS, we should emphasize that the influ¬ 
ence of carbon 2s, 2px, and 2py orbitals on the electronic 
states for the dimer-graphene system was examined. A 
comparison of the resonances with and without the inclu¬ 
sion of these orbitals showed that the contribution of C 
2s, 2px, and 2py states to the energy and intensity of the 
resonances is very small and can be ignored in our calcu¬ 
lation, indicating that the C 2pz orbitals make the main 
contribution in the emergence of the resonance states. 
Therefore, in the Hamiltonians given in Sec. II, only the 
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FIG. 4: (Color online) Calculated averaged LDOS for 
graphene in the presence of 1-5 Ca dimers at different dis¬ 
tances along the y-axis and away from center of clusters shown 
in Fig. The black dashed lines show the Fermi energy. 


C 2pz orbitals are retained in the LDOS calculations. In 
addition, we shifted all of the electronic spectra to set 
the graphene Dirac point at zero energy. In the physical 
analysis of resonances, only the energy window close to 
the graphene Dirac point has been shown. The LDOS for 
graphene in the presence of Ca dimers at different lateral 
distances from the center of each dimer-graphene cluster 
(see Fig. 1) is shown in Figs. 4(a)-4(e). Note that in 
all clusters, the central dimer is the first Ca dimer that 
is placed on the graphene monolayer and the rest of the 
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FIG. 5: (Color online) Two different arrangements of six Ca 
dimers on a graphene plane. In the text, they are referred 
to as (a) the six-dimer cluster and (b) the six-dimer ring, for 
simplicity. Note that the center-center distance between two 
adjacent dimers is 12.78 A. 
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FIG. 6: (Color online) Calculated graphene LDOS in the pres¬ 
ence of (a) a six-dimer cluster and (b) a six-dimer ring at 
different distances from the centers of the graphene clusters 
shown in Fig. along the y-axis. The black dashed lines 
show the Fermi energy. 


dimers surround it. All distances are measured from the 
center of the first Ca dimer. Figure 4(a) shows that in the 
case of one Ca dimer on the graphene layer, a weak res¬ 
onance and oscillation appears in the energy window of 
the plot; the electronic states above the Fermi energy are 
affected the most by the presence of the Ca dimer. These 
Ca-related features decrease in strength with increasing 
distance from the Ca dimer (see the LDOS in Fig. 4(a) 
at 8.48 nm). The LDOS also displays an asymmetric be¬ 
havior around the Dirac point energy that persists even 
at large distances from the Ca dimer, partly due to the 
inclusion of overlap matrix in onr theory. In this 

study, to find the Fermi energy for each dimer-graphene 
system, the energy of highest occupied molecular orbital 
(HOMO) was first determined for the central graphene 


ribbon with and without Ca dimers. Then the difference 
between the two HOMO energies was added to the Dirac 
point energy as our estimate of the Fermi energy for the 
dimer-graphene system. This procedure was repeated for 
all of the dimer-graphene clusters. The electrochemical 
potential in the graphene leads was taken to be the Dirac 
point energy renormalized to match the Fermi energy as 
determined according to the above procedure. 

In the case of the two-dimer cluster [Fig. 4(b)], the 
oscillation in the electronic LDOS above the Dirac point 
energy is stronger and has shifted downwards in energy. 
Moreover, the lowest energy resonant state is now partly 
populated with electrons, signaling that the atomic col¬ 
lapse is beginning to take place. This situation resembles 
the onset of atomic collapse in the STM experiments m, 
but corresponds to the case of four Ca dimers adsorbed 
on graphene in the experiment. 

Possible reasons why the onset of atomic collapse oc¬ 
curred for a larger number of Ca dimers in the cluster in 
the experiment of Wang et al. EOl than in the present 
theory may include the following: In the experiment [30], 
the graphene layer was on a BN substrate, whereas in 
the present work, the calculations were carried out for 
suspended graphene in vacuum (with no substrate). Di¬ 
electric screening due to the substrate in the experiment 
is expected to reduce the strength of the Coulomb po¬ 
tential well induced in the graphene by the Ca dimers, 
thus requiring a larger number of adsorbed dimers for 
the atomic collapse to take place. The larger spatial sep¬ 
aration between the Ca dimers in the experiment than in 
the present theory is also expected to result in a weaker 
Coulomb well in the experimental system, with qualita¬ 
tively similar consequences. 

As the number of Ca dimers increases from two to five 
[Fig. 4(b)-4(e)], the lowest resonant state becomes com¬ 
pletely occupied and the number of quasi-bound states 
seen in the energy window in Fig.4 increases. In fact, a 
series of resonances develops which may be related to the 
Bohr-Sommerfeld quantization as was discussed in Ref. 
|6]; however, the resolution of the STM experiment [30] 
was only able to detect the lowest state. As can be seen 
in Fig. 4, the energy of the lowest atomic collapse state 
shifts below the graphene Dirac point as more dimers are 
added to the group, while the Fermi level moves up in 
energy. This means that additional charge transfer from 
Ca atoms to graphene occurs in the process of adsorbing 
each additional Ca dimer on the graphene monolayer. 

In the experiment of Wang et al. [30| , a gaplike fea¬ 
ture in normalized differential conductance was reported 
which comes from the effects of inelastic tunneling and 
the lifetime broadening by electron-phonon and electron- 
electron interactions mi- Since we are only interested 
in studying the resonance features coming from the Ca 
dimers in suspended graphene LDOS, these effects have 
not been included in our calculations. 

To examine the spatial dependence of the resonant 
states, the electronic LDOS were computed at different 
distances from the center of each dimer-graphene cluster. 
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The results are presented in Fig. 4. Although the inten¬ 
sities of the resonances in the LDOS depend strongly on 
the distance from the center, the effect of resonances on 
the electronic spectra is still discernible even 8.48 nm 
from the center (the black curves in Fig. 4). This sug¬ 
gests that the resonant states are located mainly on the 
graphene as distinct from the localized atomic orbitals 
on Ca atoms. 

To prove this point definitively, we removed the Ca 
atomic orbitals from the our model tight-binding Hamil¬ 
tonian but retained the EPE values on the graphene that 
were obtained in our DET calculations for the systems 
consisting of the graphene and adsorbed Ca dimers. In 
this way, only the effect of charge transfer between the Ca 
atoms and graphene on the graphene electric potentials 
was included in this test calculation, while the role of Ca 
valence orbitals was ignored. Interestingly, the results 
showed a negligible change in the electronic density of 
states when compared with the results in Eig. 4, show¬ 
ing that the presence of the resonance states is mainly 
due to the charge transfer from Ca dimers to graphene 
and not due to the hybridization between Ca As and Ap 
orbitals and graphene 2pz states. 


C. Resonances in a 6-dinier cluster and ring 

We explored the dependence of the Dirac point res¬ 
onance features on the spatial arrangement of the Ca 
dimers adsorbed on the graphene by comparing the calcu¬ 
lated LDOS for two different structures, each composed 
of six Ca dimers on the graphene monolayer. The dimer 
arrangements of the two structures considered are shown 
in Eigs. 5(a) and 5(b). If we add two more dimers to 
the left and right side of the central dimer in Eig. 1(e), 
the six-dimer cluster of Eig. 5(a) is obtained. On the 
other hand, if the central dimer in Eig. 5(a) is shifted 
down along the diagonal, the 6-dimer ring shown in Eig. 
5(b) is formed. The calculated local densities of states 
for the six-dimer cluster and the 6-dimer ring are shown 
in Eigs. 6. The LDOS for the 6-dimer structures is based 
on the EPE values whose averages are shown in Eig. 3(c). 
The averaged EPE value for the six-dimer ring does not 
change significantly as a function of R for R < 6A. This 
is related to the absence of a central dimer in that struc¬ 
ture and is a completely different behavior than what is 
seen in Eig. 3(c) for the other dimer-graphene clusters. 

Eor both of the six-dimer structures the potential well 
in the range of R > 7.5k is deeper than that for clus¬ 
ters composed of 1-5 Ca dimers, indicating more quasi¬ 
bound states for structures of this type. Eigures 6(a) and 
6(b) clearly show two distinct quasi-bound states near 
the Eermi level. Interestingly, the LDOS in the two ar¬ 
rangements shows that the second atomic-collapse state 
is becoming populated. We note that only a single atomic 
collapse resonance was observed to be populated in the 
STM experiment of Wang et al. [30]. However, Wang et 
al. ED] did not report any studies of six-dimer clusters. 


Although some of the the resonance features are some¬ 
what narrower for the six-dimer ring than that for the 
six-dimer cluster in Eig. 6, it is evident that the ener¬ 
getic and spatial characteristics of the three lowest LDOS 
resonances are not affected greatly by the change in the 
arrangement of the dimers. The most significant differ¬ 
ence between the LDOS for the two geometries is in the 
resonance energies above 0.2 eV at short distances (2.21 
nm) from the center of clusters. 

These results suggest that the adsorption of more Ca 
dimers on the graphene monolayer results in the transfer 
of more electric charge from Ca atoms to the graphene 
monolayer and a deeper electrostatic potential well in 
the graphene and that, accordingly, more atomic collapse 
states will be populated. 


IV. CONCLUSION 

In summary, we have studied theoretically the 
Coulomb bound states and Dirac point resonances in¬ 
duced in a suspended graphene monolayer due to the 
adsorption of different number of Ca dimers on the 
graphene. We performed ab initio DET calculations of 
the relaxed geometry for graphene with a single adsorbed 
Ca dimer and also the Mulliken atomic charges and the 
atomic electrostatic potentials for 1-6 Ca dimers on the 
suspended graphene. The tight-binding extended Hiickel 
model, modified to include the electrostatic potentials 
obtained with DET, was then used to compute the elec¬ 
tronic structure of the system consisting of the adsorbed 
Ca atom dimers and the carbon atoms of the graphene 
monolayer. Our findings indicate that a charge trans¬ 
fer from dimers to the graphene occurs and increases as 
the number of dimers in the adsorbed cluster is increased 
from one to six. 

Our calculated electronic structures reveal that the 
first atomic-collapse state close to the Dirac point en¬ 
ergy begins to become populated with electrons when 
two Ca dimers are adsorbed on the suspended graphene. 
By contrast in the recent STM experiment [30] , the first 
atomic-collapse state was observed to begin to populate 
with electrons when four Ca dimers were present in the 
Ca cluster adsorbed on the graphene. Possible explana¬ 
tions of this difference may include dielectric screening of 
the Coulomb potential well in the graphene due to the 
BN substrate underlying the graphene monolayer in the 
experiment [30] and a shallower potential well due to the 
larger spacing between the Ca dimers in the experimental 
system. 

Our theory also predicts a second atomic-collapse state 
to begin to populate with electrons when six Ca dimers 
are adsorbed on the suspended graphene. No studies in¬ 
volving adsorbed Ca clusters with six Ca dimers were 
reported by Wang et al EO]. The experimental study [30] 
also only found evidence of a single atomic-collapse res¬ 
onance becoming populated with electrons. Therefore, 
experiments testing our prediction of a second atomic- 


collapse state being populated with electrons for larger 
numbers of Ca dimers adsorbed on suspended graphene 
would be of interest. We note that suspended graphene 
monolayer nanostructures, although without adsorbed 
Ca dimers, have already been realized experimentally and 
their transport properties have been measured. [42] 

The calculated spatial characteristics of the resonance 
features confirm that only the charge transfer from Ca 
atoms to graphene is responsible for the formation of the 
atomic-collapse resonances and that the local Ca valence 
orbitals are not involved significantly in the resonant 


states. Furthermore, it is found that the spatial arrange¬ 
ment of Ca dimers within the adsorbed dimer cluster does 
not strongly affect the formation of the atomic-collapse 
state and the location of quasi-bound states that are close 
in energy to the Dirac point, although the overall size of 
the Ca dimer cluster may be important. These findings 
can be used in future nanoscale devices where there is a 
need to confine electrical charges to a small area. 

This work was supported by NSERC, CIFAR, West- 
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